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Abstract 

We suggest a new method for the calculation of the nonlocal part of the 
effective action. It is based on resummation of perturbation series for the heat 
kernel and its functional trace at large values of the proper time parameter. We 
derive a new, essentially nonperturbative, nonlocal contribution to the effective 
Q^' action in spacetimes with dimensions d > 2. 
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1. Introduction 

X 

Effective action is among the fundamental ideas of modern quantum field theory. Cal- 
culated analytically for a given background field, it gives information about the induced 
energy-momentum tensor of quantum fields and quantum corrections to the classical 
equations of motion. The nonlocal part of the effective action should contain, for 
instance, particle creation effects. For black holes it should be able to account simulta- 
neously for vacuum polarization and asymptotic Hawking radiation. Various important 
applications of the effective action can also be found in fundamental string theory. The 
Lorentzian effective action, which we actually need, can be obtained from the Euclidian 
one, -T [</>], via analytic continuation. In turn, r[4>] can be defined by the following 
path integral 

exp(-r[#r)])= J D <pexp(-S[<p] + {<p-4)¥^\ (1.1) 

where <fi{x) is a given mean field, and the functional integration over quantum fields 
ip(x) is assumed. The general semiclassical expansion of r[(f>] begins with the one-loop 
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contribution, which is given by the gaussian path integral 

exp(-r[0(x)]) = J D(pexp(~ J dx^g<p(x) F (V, <f>(x)) <p(x) J (1.2) 

The operator F(V,(j)(x)) here determines the propagation of small field disturbances 
(p(x) on the background of 4>(x) and in bosonic case can generically be written down 
as 

F = -D + V(x), (1.3) 

where □ = V 2 = V M V u is the Laplacian in the Euclidean field theory, which becomes 
the d'Alembertian when analytically continued to the Lorentzian sector, and V(x) is 
the potential term. Note that, for some fields, the one- loop contribution is exact, for 
instance, for the scalar field without self-coupling. For a properly defined measure the 
gaussian integral (|1.2|) can be formally calculated as 



where A are the eigenvalues of the operator F corresponding to appropriately nor- 
malized eigenf unctions (f)\(x), f dx ^fg cj)\(x)(j)y(x) = 5\\>. Here the functional trace 
Tr does not depend on a particular basis in the functional space of disturbances ip 
and, therefore, in an appropriate representation it reduces to the integral over spatial 
coordinates x of the diagonal element of the operator kernel. 

The effective action (|1.4j ) is, of course, ultraviolet divergent and should be regu- 
larized, with the subsequent interpretation of explicitly isolated divergences in terms 
of infinite renormalizations of the coupling constants of the theory. These divergences 
are well understood and it is unlikely that anything new can be added here. There- 
fore, we concentrate on more interesting finite, and generally nonlocal, contributions 
to one-loop effective action. These contributions depend on infrared properties of the 
theory and contain nontrivial information about real physical effects. Analytical cal- 
culational schemes for r are usually based on the following integral representation of 
the functional trace of F 

oo 

Tr lnF = - / — Tre" sF , (1.5) 
o 

where all local divergences can be easily isolated with the aid of dimensional regular- 
ization. The kernel 

K(s\x,y) = exp(-sF)5 {d \x,y), (1.6) 
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where d is spacetime dimensionality, obviously satisfies the heat kernel equation 

j-K(s\ x, y) = -FK(s\ x, y) = (n - V{x))K(s\ x, y), (1.7) 

with the initial condition 

K(0\z,y) = 6W(z,y) (1.8) 

at s — 0. The auxiliary parameter s is usually called the proper time. Thus, calculating 
the effective action can be reduced to the solution of the Cauchy problem for K(s\ x, y). 
In fact, what we actually need is the coincidence limit of this function, since in the 
representation we used, the functional trace of the operator e~ sF corresponds to the 
integration of the diagonal elements of K over the spacetime coordinates x, so that 

r=-~Jdx(J™^-K(s\x,x) S J. (1.9) 

It is clear that the success of the calculation of mainly depends on our ability to 
find an analytical solution of the heat kernel equation and carry out the integration 
over the proper time in ( p..9|) . The integral is obviously divergent as s — > 0. As we 
have already mentioned above, this divergence can be easily isolated and interpreted in 
terms of the local ultraviolet properties of the theory. On the other hand, the behavior 
of the integral at infinity, s — > oo, determines infrared properties of the theory and 
carries the physical information, e.g., particle creation. If the field has a big positive 
mass the proper time integral is convergent at s — > oo. However, in the case of massless 
fields, the situation is much less trivial. The infrared convergence here depends on the 
approximation scheme used to calculate K(s\x,y). Then, it is often even unclear to 
what extent the obtained effective action reflects the physical properties of the theory 
rather than the features of the approximation scheme used. 

Below, we discuss the known calculational techniques, namely, the local Schwinger- 
DeWitt expansion fl], ||, nonlocal covariant perturbation theory 0, £|, ^[ and the modi- 
fied gradient expansion || and point out why all of them fail when applied to interesting 
physical problems. Instead of them, we suggest a new method based on resummation 
of perturbation series and calculate new, essentially nonperturbative terms in the ef- 
fective action. This method becomes indispensable in low- dimensional models (d < 2) 
where all previously known techniques are inapplicable. In this paper, we demonstrate 
how our method works in flat space of dimension d > 2, while the generalization to the 
curved space and low- dimensional case will be considered in 0. 

One of the main results of this paper is an exact (nonperturbative in V) late- 
time asymptotics for the heat kernel which in a spacetime of dimension d > 2 for the 
coincidence limit takes the following form 

1 



K(s\x,x) = j—^ J - 2 i [ l + w - v V(x)) , s -oo. (1.10) 
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This asymptotics can also be found when the arguments of K(s\x,y) are different - 
see Sect. 4 for details. Another important result is the nonperturbative in potential 
expression 

TrK(8) = j^Jdx{l-s(v(x)+V^V(xj)} s ^ oo (1.11) 

for the functional trace of i^(s)Q. 

To avoid excessive use of integration signs we employ here and throughout the paper 
the following shorthand notation 

J(x)= [ dyG{x,y)J(y), (1.12) 



n-v 



where G(x, y) is the Green's function of the operator — F = □ — V with zero boundary 
conditions at spacetime infinity, that is 

(D-V)G(x,y) = 5^(x,y), G(z,y)-0, |x| - oo, (1.13) 

and J{x) can be any function of various field quantities like powers of potential, its 
derivatives, etc. We always presume that the spacetime has the positive definite (Eu- 
clidean) signature, so that the Laplacian □ is negative definite assuming zero boundary 
conditions at infinity. Moreover we consider only non-negative potentials V (x) > 0, 
so that the whole operator F = — □ + V is positive definite. Therefore, the Green's 
function ( J1.13| ) is uniquely defined and guarantees that the nonlocal expression ( |1.12j ) 
makes sense for d > 2. 

As we shall see, the asymptotics ( l.lOf) and Ql.ll ) are the corner stone of the 



technique we develop for the calculation of nonlocal contributions to the effective action. 
In particular, they lead to essentially nonperturbative terms which can be explicitly 
calculated for two broad classes of potentials with compact support, namely, for those 
ones which are, respectively, very small or very big in units of the inverse size of their 
support. For small potentials we get the terms which replacethe conventional Coleman- 
Weinberg contribution to the effective action. In four dimensions, for instance, these 
terms read 



1.14) 



where the mass parameter /x 2 reflects the usual ultraviolet renormalization ambigu- 
ity. On the contrary, in the case of big potentials the Coleman- Weinberg action is 

-"^Note that this expression for Tr K(s) cannot be obtained directly by integrating the asymptotics 
( 1.10| ) over the whole spacetime because for a given s this asymptotics fails at \x\ 2 > s. Its derivation 



is given in Sect. 4 and Appendix B, where we show that the expressions ( 1.10 ) and ( 1.11 ) are in 
complete agreement with each other. 
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supplemented by the nonlocal term of the form 

Ar = 6^ / **( v + v nh v Y- < L15 > 

\x\<R 

Here R is the size of the compact support of V(x), that is, V(x) = at \x\ > R, and 
(...) denotes the spacetime averaging of the corresponding quantity over this compact 
domain. The obtained expressions are both nonlocal and non-analytic in the potential 
V(x). 

The paper is organized as follows. In Sect. 2 we consider the known approximation 
schemes and discuss their applicability in the infrared region. Sect. 3 is devoted to 
the nonlocal and nonlinear resummation of the Schwinger-DeWitt perturbation series, 
corresponding to the so-called connected graph expansion of the heat kernel. In Sect. 4, 
as an extention of this resummation, we derive the asymptotics (|1.10| ) and discuss 



its relation to the functional trace of the heat kernel. The nonperturbative, 

nonlocal contributions to the effective action are obtained in Sect. 5 with the aid of 
the new technique based on a piecewise smooth approximation for the heat kernel. 
In two appendices we give details of the resummation technique and derive from the 
covariant perturbation theory the asymptotics of the heat kernel trace ( |1 . 1 1| ) up to the 
first subleading order in 1/s inclusive. 



2. Approximation schemes and infrared properties 
of the effective action 

In flat space, which we consider in this paper, the solution of the heat kernel equation 
can be easily found if the potential vanishes. For an arbitrary spatially dependent 
potential the analytical expressions are, of course, available only in certain approxima- 
tions. In the general case, it is convenient to factorize the "zero potential" part of the 
solution explicitly and use the following ansatz for K(s\ x,y) 



K(s\x,y) 



1 



(Aits 



A 12 



exp 



\x - y\ 
4s 



Q(s\ x,y), 



(2.1) 



where the factor singular in s guarantees that the initial condition ( |1.8| ) is satisfied, 
provided that Q is analytic in s at s = and Q(0 \ x, y) = 1. If V = 0, then = 1 and 
hence all nontrivial information about the potential is encoded in the deviation of Q 
from unity. 

The most well-known approximation used for the calculation of K(s \ x, y) is the so- 
called local Schwinger-DeWitt expansion, where Q is written down as a series in growing 
powers of the proper time s. This expansion is a very powerful tool for revealing local 
ultraviolet properties of the theory. However, when applied in the infrared region, it 
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gives a finite result only for massive fields. If the potential V(x) has a large positive 
constant part, that is, 



V(x) = m 2 + v(x), 



(2.2) 



where m? is the squared mass of the field, then the function Q(s \ x, y) contains an overall 

2 

exponential factor e~ sm and the independent of mass part is expanded in powers of s 



Q(s\ x,y) = e sm2 a n (x, y) s 1 



(2.3) 



n=0 



Here a n (x, y) are the two-point Schwinger-DeWitt coefficients, whoes coincidence limits 
(x — > y) are explicitly calculable in general field theories, including gravity. Substitut- 
ing ( |2.3|) in (|2.1|) and then the obtained expression at x = y in ( |1.9|) one gets: 



r 



i 



2 (An) 



d/2 



/*£ (I 



s n~d/2-l e - 



ds ) a n (x, x). 



(2.4) 



It is important that the exponent e~ sm is not expanded here in powers of s. There- 
fore, in the proper time integral it provides a cutoff at the upper limit, so that the pow- 
ers of s in this expansion get effectively replaced by powers of 1/m 2 . The first (d/2 + 1) 
integrals in ( |2.4| ) diverge at s — > and should be regularized. To do that, we apply 
the dimensional regularization method. Namely, by replacing the dimensionality d by 
2a;, we calculate the integrals in the domain of their convergence and then analytically 
continue the result to uj — > d/2. In spaces with even number of dimensions, which we 
mainly consider in what follows, this gives rise to the contribution -Tdiv, log containing 
the pole at u = d/2 and the term logarithmic in m 2 



i 



div,log 



2(4tt 



d/2 

dx yj - 

n=0 



—m) ■ 



(f-n)! 



n + 1 +ln 



m 

47T/i 2 



(2.5) 



where uo — > d/2. The pole corresponds to an infinite ultraviolet renormalization of 
the terms proportional to ao,---,ctd/2 m the original Lagrangian. Other terms in the 
expansion (|2.4|) are finite and give the infrared contribution to the total action 



r — -r d i V . 



in 



m 2 \ 
4^J 



d/2 



dx 

n=d/2+l 



T(n-d/2) 



[m 



2\n 



(2.6) 



The Schwinger-DeWitt coefficients the homogeneous polynomials of di- 

mensionality 2n in the units of inverse length, which are built of v{x) and its multiple 
derivatives . Therefore, on dimensional grounds, they can be symbolically written down 

as 

a n (x,x) ~ v k (x)(V l v j )(x), 
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where i denotes an overall number of derivatives acting in all possible ways on j factors 
of v(x), and k powers of v(x) stay undifferentiated. The positive integers (k,j,i) are 
related to n as 2(k + j) + % = In. It is clear that the infinite series in ( |2.6|) represents 
the expansion in growing powers of the following dimensionless quantities 

v (x) , V l v(x) 

-V<1, — (2.7) 

which obviously should be much smaller than unity. Only in this case the first few 
terms in the asymptotic series Q2.6Q are reliable. 

Thus, the Schwinger-DeWitt expansion is applicable only in theories with small and 
slowly varying fields as compared to a big mass parameter. This expansion contains 
only local terms. This is not surprising because all nonlocal effects, e.g., particle 
creation, are very small for heavy particles in a weak external field and cannot be 
handled by this method. The Schwinger-DeWitt technique can be easily extended to 
curved spacetime and to theories with covariant derivatives built with respect to an 
arbitrary fibre-bundle connection. In this case, the perturbation potential v(x) will also 
depend on the spacetime curvature tensor and fibre-bundle curvatures (commutator 
of covariant derivatives). The smallness of fields and their derivatives includes the 
requirement of the smallness of these curvatures and their derivatives as well. Despite 
its universality, the Schwinger-DeWitt expansion becomes inefficient when the ratios 
in ( |2.7| ) become of the order of unity and completely fails for massless fields. In the last 
case all integrals over the proper time integral are infrared divergent. This divergence 
has, of course, no physical meaning and is an artifact of the approximation technique 
used. 

There are two known ways to proceed with massless fields. One possibility is the 
resummation of all terms which contain the undifferentiated potential V(x) in the local 
Schwinger-DeWitt series (|2.3|). They are summed up to form an exponent similar to 



Q(s\ x, x) = e~ sV{x) ^(xi x) s n . (2.8) 

71=0 

This method was suggested in |J, where a regular technique for the calculation of the 
modified Schwinger-DeWitt coefficients a n (x,y) was also presented. The proper time 
integral in ( |1.9|) has now an infrared cutoff at s ~ 1/V(x) and in this case the effective 



action is similar to ( |2.5| ) - (|2.6|) , where m is replaced by V(x) and a n (x, x) by a n (x, x). 



It is convenient to write this action as a sum of three terms 

r = r div + r cw + r fin , (2.9) 



where the divergent part is equal to 



, , d 



Q"n \X; X J . 

(2.10) 
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The pole part of this action coincides with that of (|2.5|) if we take m 2 — > limit of 
(|2.5|). Actually, in this case, only the term proportional to survives in -Tdiv.iog and 
by virtue of the relation between twiddled and untwiddled coefficients, namely, 



d/2 



d/2-n 



n=0 



the pole parts of (|2.5|) and ( p.lOj ) are the same. The terms proportional to T'(d/2 — n+1) 
perform finite renormalization of the local terms y d l 2 - n a n . The logarithmic terms of 
( |2.5|) are replaced in the modified action (|2.9| ) by 



1 r ^(_yW2-n V 1 f V 

rcw = ofA \dn / dx 2^ 777^ V\ ln — = ofA w/ 2 / ax In — a d/2 . 

2(47r) d / 2 J f^(d/2-n)! /i 2 2{4n) d / 2 J ^ 



n=0 



This is nothing but the space-time integral of the Coleman- Weinberg effective potential. 
For instance, in four dimensions the leading term is the original Coleman- Weinberg 
effective potential, V 2 \ii(V/ij, 2 )/6A7t 2 , while the rest represents corrections due to the 
derivatives of V(x). Similarly to (|2.6| ), the finite part Jfi n is an infinite series 

J V 7 n=d/2+l V ' 

The modified Schwinger-DeWitt coefficients do not contain the undifferentiated po- 
tential and the typical structure of the terms entering a n (x, x) is V m ^ J (x), where 
m + 2j = 2n. Every V here should be differentiated at least once and therefore m > j. 
Thus the coefficients be symbolically written down as 

[2n/3] 

fi n (x,x)~ ^ 2n ~ 2:i V\ (2.14) 

3=1 

where the upper value of j is the integer part of 2n/3. 

This perturbation theory is efficient as long as the potential is slowly varying or 
bounded from below by a large positive constant, so that 

™ «!,««! ,2,5) 

V 2 (x) V 6 (x) 

The case of bounded potentials reproduces the original Schwinger-DeWitt expansion for 
nonvanishing mass. Therefore, let us consider the potentials which vanish at spacetime 
infinity (|x| — > oo). Namely, we assume the case of power-like falloff 
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for some positive p. For such a potential terms of the perturbation series (|2.13 ) behave 

as 



[2n/3] 

y- n ~ E \4 p - 2)(n - j) (2-17) 
j'=i 

and thus decrease with n only if p < 2. For p > 2, the modified gradient expansion 
completely breaks down. It makes sense only for slowly decreasing potentials of the 
form ( |2.16| ) with p < 2. In this case the potential V(x) is not integrable over the 
whole spacetime (J dxV(x) = oo) and moreover even the operation (l/0)V(x) is not 
well defined^. Therefore the above restriction is too strong to account for many in- 
teresting physical problems. In addition, similarly to ( [2.6|) , the asymptotic expansion 



(|2.13|) is completely local. It does not allow us to capture nonlocal effects, which are 
exponentially small for potentials satisfying (|2.15|) . 

The way to overcome this difficulty and to take into account nonlocal effects was 
suggested in the form of covariant perturbation theory (CPT) |], In this theory 
the full potential V(x) is treated as a perturbation and the solution of the heat equa- 
tion is found as a series in its powers. From the viewpoint of the Schwinger-DeWitt 
expansion it corresponds to an infinite resummation of all terms with a given power of 
the potential and arbitrary number of derivatives. The result reads as 

/oo 
dxK(s\x,x) = J2 TrK n(s), (2-18) 
n=0 



where 



Ti K n (s) = J dx 1 dx 2 ...dx n F n (s\xi,x 2 , ...x n )V(xi)V(x 2 )---V(x n ), (2.19) 

and the nonlocal form factors F n (s\ x\,x 2 , ...£«) were explicitly obtained in ]3|, fl], |5|. It 
was shown that at s — > oo the terms in this expansion behave as 

Tr K n (a) = O (p^rr), n>l, (2.20) 

and, therefore in spacetime dimension d > 3 the integrals in (|1.9p are infrared conver- 
gent 



s \s 



In one and two dimensions this expansion for r does not exist except for the special 
case of the massless theory in curved two-dimensional spacetime, when it reproduces 



2 For the convergence of the integral in (1/n) V the potential V(x) should fall off at least as l/|a:| 3 
in any spacetime dimension B. 
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the so-called Polyakov action [|, |[ UfJ CPT should always be applicable whenever 
d > 3 and the potential V is sufficiently small[], so that its effective action explicitly 
features analyticity in the potential at V = 0. Therefore, its serious disadvantage is 
that this theory does not allow one to overstep the limits of perturbation scheme and, 
in particular, discover non-analytic structures in the action if they exist. 

All this implies the necessity of a new approximation technique that would allow us 
to overcome disadvantages of the above methods. In the rest of this paper we develop 
such a technique, involving further resummation of the perturbation series. We develop 
an infrared improved perturbation theory for the heat kernel and reveal new nonlocal 
and non-analytical structures in the effective action. 



3. Resummation of proper time series 

We use the exponential ansatz for the function fi(s| x,y) defined by Eq. (|2.1|) : 

Q(s\ x, y) = exp — W(s\x,y) 

Our goal is to develop an approximation technique for W similar to CPT, which is an 
alternative to the expansion in s. By virtue of (|1.7|) and ( |1.8| ) the function W(s\ x,y) 
satisfies the equation 



dw {x - y y 



ds s 
with the initial condition 

W{s = 0\x,y) = 0. 



v^w -uw — V— (yw) 2 , 



(3.1) 



(3.2) 



This equation is nonlinear in W and we solve it by iteration, considering (VII 7 ) 2 as a 
perturbation. For this purpose it is convenient to rewrite ( |3.1| ) - ( |3.2j ) as an integral 



equation. In Appendix A, it is shown that this integral equation takes the following 
form 



W(s\ x, y) 



dae sa{1 ~ a) ° (V(x) - (VW(sa\x,y)) 2 ) 



(3.3) 



where the operator □ acts on the argument x = x(a\x, y) = ax+(l— a)y. The equation 
(|3.3|) can be easily solved if (VW) 2 V. As we will see later, this condition is satisfied 

3 This action can be obtained by integrating the conformal anomaly |)| ||. 

4 The conditions of the smallncss of the potential are exactly opposite to those of (2.15), e.g., 



V 2 /V 2 V <C 1. However, this is true only as a rather rough estimate, because CPT is a nonlocal per- 
turbation theory and its actual smallness "parameters" are some nonlocal functionals of the potential. 
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for a broad class of potentials V. The lowest order approximation for W is obtained 
by omitting (Vlf) 2 in eq. (|OD 



Wo(s\ x, y) = s I da e^ 1 -) 5 V(x) \ s=ax+{1 . a)y ■ (3.4) 



o 



This is a linear but essentially nonlocal functional of the potential. Further terms of the 
perturbation theory, W n = O [(VWo) n ], can be graphically represented by connected 
tree graphs with two derivatives in the vertices, internal lines associated with the 
nonlocal operator 

f(-sn) = [ dae sa ^ D , (3.5) 
Jo 

and external lines given by ( |3.4Tj f|. Note that this connected graph structure arises in 
the exponential and when expanded gives rise to the disconnected graphs. In context 
of the heat kernel expansion this property was observed in |TIJ . Resummation of the 
perturbation series in V explicitly features exponentiation of the quantities contain- 
ing only connected graphs. Here we showed how this exponentiated set of connected 
graphs directly arises from the solution of the simple nonlinear equation (|3.1|) . The 
"propagator" (|3.5| ) was worked out within the covariant perturbation theory in |3], |], [5] 
and was also obtained in || by direct summation of gradient series. 

At this stage the efficiency of the connected graph expansion is not yet obvious. 
Crudely, it runs in powers of the dimensionless quantity (s/(— sD)V) 2 V(x) which, at 
least naively, should be small for slowly varying or/and small potentials. Apart from 
this, infrared properties of the effective action strongly depend on the lowest order 
approximation for W fl3.4|) . The effective action involves only diagonal elements of the 
two-point function W(s\x,y), which look much simpler than ( |3.4j) 

W (s\x,x) = s I dae sa{1 ~ a)n V(x). (3.6) 
Jo 

Note that, at small s, the function Wo can be expanded as Wo = sV + O (s 2 ) . The only 
term with undifferentiated potential entering Wq is linear in s, while all other terms 
contain derivatives of the potential V. The same is also true for the exact W, which 
differs from Wo by higher powers of the differentiated potential. This completely agrees 
with the modified gradient expansion discussed in Sect. 2. However, the expression (|3.6j ) 
directly involves the nonlocal operator and its late time behavior is very different from 
that naively expected in the modified gradient expansion and in CPT. 

To show that, let us first find the coordinate representation of the operator (|3.5|) , 
that is, the kernel f(—sO)5^(x,y). Using a well-known result for the exponentiated 

5 This graphical interpretation should not be taken too literally because integration over a- 
parameter(s) involves also the argument x = x(a\x,y) of the potential. 
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□-operator 



sa(l— a)CI 



5 {d) (x,y) 



(47rsa(l - a)) d / 2 



exp 



a) 



(3.7) 



one can write 

dae sa(1 ~ a)D 5^(x,y) 

p -\x-y\ 2 /2s 



(47Ts) rf / 2 



d(3 



d-2 



fjd/2 



exp 



x 



As 



(3.* 



where we have changed the integration variable, at = + < /3 < oo. For d > 3, 
the integral can be easily calculated and the result is expressed as a sum of McDonald 
functions of the argument \x — y\ 2 /2s 



f{-sn)8W(x,y) 



9 — \x— y | 2 /2s d— 1 



k=l 



where C%_ 2 are the binomial coefficients. For very large s, the argument of K k 



(3.9) 



dlZ is 



small. Using the asymptotics: ^(2) — T(\v\){2/ z)^ /2, z 
s the form factor is dominated by the following term 



f(-sn)S^(x,y) 



1 r(d/2-i; 



2svr d / 2 |x 
Taking into account the fact that 



Id— 2 



+ 



0, we find that at large 



(3.10) 



□ * W (s,v) = 
we finally obtain: 



r(d/2 - 1) 



4vr d / 2 



y 



d-2 ' 



W (s| x, x) = sf(-sn)V = V(x) + O 



□ 



(3.11) 



(3.12) 



This behavior agrees with the formal asymptotics found by the Laplace method in |4[]. 
Therefore, at s — > 00, the function Wo approaches a constant. The nonlocal functional 



00. 



(|3.12|) is well-defined for d > 3 only if the potential V vanishes fast enough at 

For splitted arguments the asymptotics of Wo(s| x, y) is more intricate. In this case 
the form factor ( |3.5| ) no longer arises as a whole because the integration parameter 
a appears in eq.( |3.4| ) also in the argument x = ax + (1 — a)y of the potential V{x). 
Applying the Laplace method, one can show that the integral ( |3.4| ) is dominated by 
the contribution of the end points, a = and a — 1. These contributions are different 
because x(a = 0) = y and x(a = 1) = x, whence 



W (s\x,y) 



□ 



V{x) 



□ 



V(y) + O 



(3.13) 
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Substituting this expression in ( p.3|) and solving the integral equation by iterations one 
can find the late time asymptotics W^x, y) for the exact W(s\ x, y) 

W(s\ x, y) = W^x, y) + (-) (3.14) 



as a nonlocal gradient series 

W 00 (x 1 y) = ~V(x) + ^(^VV(x)\ h... + (x^y). (3.15) 



□ v ' □ V □ 



It is remarkable, however, that this series can be "summed up" and the nonlocal 
expression for Woo(x,y) can be found exactly in terms of the Green's function of the 
original operator F — — □ + V. 



4. Late time asymptotics of the heat kernel 

By substituting the late time ansatz ( |3.14j ) in equation (|3.1|), it is easy to see that the 
first two terms in its left hand side vanish at s — > oo, while the rest reduce to the 
equation for (x,y) 

(□ - V) e - w °° {x ' y) = 0. (4.1) 

Despite the positivity of the operator — □ + V with Dirichlet boundary conditions at 
\x\ — > oo, this equation admits non-trivial solutions. In fact, e~ Woo( - x ' y ^ does not have 
to go to zero at infinity^]. In view of the iterative solution ( |3.15| ) it should tend to some 
unknown function of y 

The equation ( |4.1| ) with this boundary condition is then solved by 

e -W x {x,y) =C (y)$(x), (4.3) 

if the new function satisfies the equation 

(□ - V) = 0, (4.4) 

and $(a;) — > 1 at \x\ — > oo. The solution of this problem for $(a;) is uniquely determined 
in terms of the Green's function ( |1.13| ) 

$( x ) = l + —L-V(x). (4.5) 



6 The boundary condition K(s\x,y) — > at \x\ — » oo is enforced by the gaussian factor in ( ^.l| ), 
even for nonvanishing finite f2(s| x, y) = cxp(— W(s\ x, y). 
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This function of x is a nonlocal and essentially nonlinear functional of the potential, 
and it will play a very important role in what follows. 

The heat kernel is symmetric in its arguments x, y and, therefore, the unknown 
function C(y) should coincide with Thus, finally we obtain the following exact 

late time asymptotics 

e -W»(*.w) = $(a;)$(y). (4.6) 

Expanding Woo(x, y) in powers of the potential V one gets 

Woo{x, y)=- In $(x) - In $(y) 

1 ,„ , 1„ 1 ,„ , 1/1 N 2 



=- 5 V(ar) - -^-V^(x) + - I -V(x) ) + ... + (x <- y). (4.7) 



The first term here is in agreement with the perturbative asymptotics fl3.15| ). However 



beyond that, the iterative solution (|3 .15 ) seems to be in contradiction with (4.7|). The 



series ( J3.15|j runs in powers of the differentiated potential, while the expansion ( |4.7| ) 
contains only undifferentiated potential. However, integration by parts in the second 
term of (|3.15[) exactly reproduces the second and third terms of ( |4.7|) . Thus, both 
expansions are equivalent, but the first one reveals more explicitly the gradient of the 
potential as a small parameter, while in (|4.7Q the smallness is a result of non-trivial 
cancellations between different terms. 

Finally we write down the exact late time asymptotics for the heat kernel, advocated 
in Introduction, 

K{s\x,y)= , 4 w /2 $(g)$(y), s^oo. (4.8) 

Its heuristic interpretation is rather transparent. The heat kernel can be decomposed 
in the series 

K(s\ x,y) = J2 e- XS( ^x(x) <S> x (y), (4.9) 

A 

where A and $a are, respectively, the eigenvalues and eigenf unctions of the operator 
F — — □ + V(x). Since F is a positive-semidefinite operator, only the lowest eigen- 
mode with A = survives in this expression in the limit s —>■ oo. The appropriate 



eigenfunction satisfies the equation F§x=o — 0, which coincides with (|4.4| ) and there- 
fore &\=o(x) = $(x). The spectrum of the operator is continuous and the eigenmodes 
are not square integrable (<&o(x) — > 1 at |x| — > oo). This is why the integral over the 
spectrum, denoted above by J^ A , yields within the steepest decent approximation the 
power-like asymptotics, l/s d ^ 2 , rather than the exponential one. Of course, these ar- 
guments are not very rigorous. The zero mode $(x) with unit boundary condition at 
infinity, does not even belong to the continuous spectrum of modes normalized to the 
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delta function. Nevertheless, as we have seen, this particular mode gives the leading 
contribution to the late time asymptotics of the heat kernel. 

If we want to calculate the contribution to the effective action due to the late 
time behavior of the heat kernel we need its functional trace - spacetime integral of 
the coincidence limit K(s\x,x). Unfortunately, the expression ( |4.8| ) cannot be used 
directly to calculate Tr K(s) for a given s. Point is that this asymptotic expression 
taken at a fixed large s is applicable only for \x\ 2 < s and fails at \x\ 2 3> sQ. At the 
same time, when calculating the trace we have to integrate over the whole spacetime 
up to \x\ — ► oo and therefore need the heat kernel behavior for \x\ 2 ^> s. The attempt 
to disregard this subtlety and integrate the coincidence limit of ( Q|) over x results 
in a poorly defined quantity - the spacetime integral strongly divergent at infinity. 
Nevertheless, one can use the expression (|4.8|) to find Tr K(s) with the aid of the 
following somewhat subtler procedure. 

First let us write the variational relation 

S-ZKM = = J (4.10) 

where, of course, K(s) = exp[s(D — V)]. Then it follows that 

-^ = -^(.1.,.). (4.11) 

Substituting K(s\x,x) from ( |4.8| ) in the right hand side of this relation we obtain the 
following functional differential equation 

5V(x) (47rs)"/2 1 >' 1 } 

This equation satisfies the intergrability condition, because the variational derivative 

^ = 2W(lf)^ll) (4-13) 

is symmetric in x and y. Therefore ( |4.12| ) can be integrated to determine TtK(s). The 
solution subject to the obvious boundary conditions at V = reads 

Tr K(s) = jJ^p Jdx(l- sV$) , s^oo. (4.14) 



One can easily check that this expression satisfies the equation (|4.12 ) 



It is remarkable that in the covariant perturbation theory the leading and next 
subleading (in s) terms of the heat kernel trace can be explicitly calculated for s —>■ oo 
in every order of expansion in powers of V. The corresponding infinite series can be 



7 This follows from the derivation of (4.8) above, which is based on discarding the second term of 



eq.(3.1) linearly growing in (x — y)/s. 
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explicitly summed up to yield essentially nonlocal and nonlinear in V(x) expression 
for TrK(s). This is done in Appendix B to the first subleading order inclusive. The 
following very simple and concise result reads as 

Tr K(s) = J dx ji _ sV ® _ av^-l^V^ + O } (4.15) 

in terms of the function and its derivatives. As we see, it exactly reproduces the 
leading order term of (|4.14j) , 0(s/s d ^ 2 ), and also gives a nontrivial 0(l/s d//2 )-correction. 



Below we use this asymptotics for obtaining new types of the nonlocal effective action. 



5. Effective action 

The functional trace of the heat kernel is everything we need for the calculation of the 
effective action. Unfortunately, only its asymptotics are known. Namely, at small s 
one can use the modified gradient expansion ( |2.1| ), ( |2.8| ) and at large s - the nonlocal 
and nonlinear expression ( f4.15| ). The goal of this section is to unify both of these 



approximations to get the expression for the effective action which would incorporate 
both the ultraviolet and infrared properties of the theory. The calculation will be 
explicitly done in four dimensional case. The generalization to other dimensions d > 2 
is straightforward. 

The key idea is to replace Tri^(s) in (|1.9|) by some approximate function Tr K(s) 
such that the integral over s 

r = -\ J^TrK(s) (5.1) 

becomes explicitly calculable. The difference TrK(s) — K(s) can then be treated as 
a perturbation. Certainly, the efficiency of this procedure very much depends on the 
successful choice of K(s). Here we exploit the simplest possibility - namely, let us 
take two simple functions Tr^<(s) and Tr^>(s), which coincide with the leading 
asymptotics of Tr K(s) at s — > and s — > oo and use them to approximate Tr K (s) 
respectively at < s < and s* < s < oo. In turn, the value of will be determined 
from the requirement that these two functions match at s*. This will guarantee the 
stationarity of t with respect to the choice of s*, dF /ds* = 0. We will discuss the 
justification of this procedure a little later, while now let us proceed with the calculation 
off. 

At small s we use the lowest order term of the modified gradient expansion 



TrfQ(s) = — — -- / dx exp(— Vs) , s < (5.2) 
(4tts) z J 

and disregard all terms containing derivatives of the potential V. The rest of the asymp- 
totic series fl2.8|) , containing d n , will be treated by perturbations. Correspondingly, at 
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large s > s* we use the late time asymptotics (|4.14j) 
1 



TiKJs) 



dx (1 - sV<&), s > s 4 



(4vrs) 2 

The requirement of stationarity of t with respect to leads to the equation 



(5.3) 



j dx exp (-Vs#) = y (1 - s*V§), 



(5.4) 



which determines the value of s* as some nontrivial functional of the potential, s* = 
s*[V(x)]. Unfortunately this functional is not calculable explicitly in general, but 
nevertheless, as we will see below, it can be obtained for two rather broad classes of 
potentials. The action ( |5.1|) can be written down as a sum of two contributions 



r = r < + r > 



oo oo 

\ j d ^ K< {s)- 1 - j ^(Tri^-Tr #<(*)). 

s* 



(5.5) 



The first integral here has already been calculated and is given by the sum of the 
expressions (|2.10 ) and ( |2.12| ) with a = 1 and a n — 0, n > 1, which in our particular 
case give rise to 



r< — -Tdiv + -few 
1 



64vr 2 



dx 



— + 2C-3-ln4vr ] V 2 + V 2 \n^ 

2 — uj J n 



u^2, (5.6) 



where C = 0, 577... is the Euler's constant. The first term here is responsible for the 
renormalization of the original action and the second one is just the Coleman- Weinberg 
potential. The second integral in ( |5.5| ) can also be calculated exactly. Integrating by 
parts and taking into account Q5.4p , we obtain 



oo 

1 f ds 



2 j s 



64tt 2 



dx 



V® Ve~ s ' v 



+ V 2 T(0,s*V) 



(5.7) 



where T(0,x) is an incomplete gamma function, T(0,x) = dtt 1 e *, with the fol- 
lowing asymptotics 



r(o,x) 



In 



x 



X 



C, x < 1, 

X > 1. 



(5.8) 



Further steps strongly depend on the class of potentials, for which the consistency 
of the piecewise approximation (|5.2j ) - (|5.3|) should be carefully analyzed. 
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5.1. Small potential 



s I \ dxVfi— — (5.10) 



The approximation ( |5.2| ) - ( |5.3| ) is efficient only if the ranges of validity of two asymp- 
totic expansions (respectively for small and big s) overlap with each other and the 
point s* belongs to this overlap. In this case the corrections due to the deviation of 
Tr K(s) from the exact Tr K(s) are uniformly bounded everywhere and one can expect 
that ( [5.1D would give a good zeroth-order approximation to an exact result. Below 
we show that this necessary requirement can be satisfied at least for two rather wide 
classes of potentials V(x). 

The modified gradient expansion is well applicable in the overlap range of the 
parameter s if 

sVW < V, (5.9) 

(cf. Eq.( |2.15| ) with s replaced by effective cutoff s = 1/V) and the applicability of the 
large s expansion in the same domain reads as 

1 

V-D 

which means that the subleading term (quadratic in V M $) of the late time expansion 
(|4.15 ) is much smaller than the leading second term. 

To implement these requirements, let us make some simplifying assumptions. In- 
stead of the power law falloff, assume that V(x) has a compact support of finite size 
R 

V(x) = 0, |ar| > R. (5.11) 

Let us also assume that the characteristic magnitude of the potential inside its support 
is given by Vq. Then the estimate for the derivatives is obvious 

VW ~ || (5.12) 

and ( |5.9|) reads as sV /R 2 <C V , that is 

s < R 2 . (5.13) 

To find out what does the criterion ( |5.10| ) mean let us make a simplifying assumption, 
namely, that the potential V is small. In this case it can be disregarded in the Green's 
functions and !./(□ — V) can be replaced by I/O. Therefore the following estimates 
hold 

^V V(X) ~ Ln dV ¥^ V(v) ~ ~ V ° R2 < 

dxV<S> ~ V R d , 

dx V M $ - \ D V^$ ~ V^R d+4 . (5.14) 
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Roughly, every Green's function gives the factor R 2 , every derivative -1/R, integration 
gives the volume of compact support R d , etc. Applying these estimates to eq. ( p.lOj) 
we get sV R d > 1/ 2 i? d+4 , whence 

s > V R 4 . (5.15) 



Combining this with (|5.13|) one gets the following range of overlap of our asymptotic 
expansions 

R 2 > s > V R\ (5.16) 

It immediately follows from here that this overlap domain is not empty only if 

V R 2 < 1. (5.17) 

Moreover, the assumption of disregarding the potential in the Green's function is also 
justified in this case since V ~ Vq <C 1/R 2 ~ □. In other words this bound means that 
the potential is small in units of the inverse size of its compact support. 

Now let us check whether s* introduced above belongs to the overlap domain ( 5.16|) . 



Note that if it is really so, then s*V in Eq. ( |5.4j ) is much smaller than unity because 
in the overlap range one has sV ~ sVq <C R 2 Vq <C 1. Hence the exponent in the left 
hand side of ( |5.4j ) can be expanded in powers of s*V, and the resulting equation for s* 
becomes^ 

J dx (l- s*V + ^-V 2 + O ((s*V) 3 )^j = J dx (1 - s*V$) (5.18) 

Its solutions has the following form: 

JdxVjl-®) _ JdxVy^V 
S *^ 2 JdxV 2 ~ 2 JdxV 2 • (5 ' 19) 



Taking into account the estimates (|5.14 ) we see that the point ~ R 2 belongs to the 



upper edge of the interval ( ^.16| ). Late time expansions is fairly well satisfied here, 
but the small s expansion is on the verge of breakdown. At this level of generality it 
is hard to overstep the uncertainty of this estimate. There is a hope that numerical 
coefficients in more precise estimates (with concrete potentials) can be large enough 
to shift s* to the interior of the interval ( |5.16| ) and, thus, make our approximation 
completely reliable. 

Bearing in mind all these reservations let us proceed with the calculation of the 



effective action. Using the small x asymptotics (|5.8|) in the expression (|5.7| ) we get 

(5.20) 



s.V 



8 Note that the quadratic term should be retained in the expansion of e s " v if we want to get a 
nontrivial solution for s*. 
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It is interesting to note that in the whole action r = T < + r> the Coleman- Weinberg 
term disappears and the final answer reads 

^^\-2 L ^ + C - 2 -^)S iixV ' 
+ 64^ / j dxV2 )-e^J **V«]n(/ M^V). (5.21) 



The first term here differs from i~div i n (|5.6| ) by a finite renormalization of the local 
V 2 -term, while the two other terms have entirely new nonlinear and nonlocal structure 
advocated in the Introduction. The ultraviolet renormalization mass parameter fi 2 
makes the argument of the second logarithm dimensionless - it plays the same role 
as for the Coleman- Weinberg potential, but now it enters the new essentially nonlocal 
structure. 

It is natural that the original Coleman- Weinberg term for the case of small poten- 
tials (5.17) gets replaced by the other qualitatively new nonlocal structure. Potentials 
which are small in units of the inverse size of their support are qualitatively very 
different from nearly constant potentials for which the Coleman- Weiberg action was 
originally derived. In the case of small potentials spacetime gradients dominate over 
their magnitude and, therefore, one should not expect that the Colman- Weiberg term 
would survive the inclusion of nonlocal structures. 



5.2. Big potential 

Remarkably, the case of the small potentials ( 5.1 7|) is not the only one when one can 
find a non-empty domain of overlap where both asymptotics for Tr K(s) are applicable. 
Namely, the opposite case of big potentials (in units of the inverse size of their support) 

V R 2 > 1, (5.22) 

is equally good. The key observation here is that in this case the kernel of the Green's 
function !./(□ — V) can be replaced within the compact support of by — 1/V (□ ~ 
1/R 2 <C V Q ~ V) and correspondingly 

' V(x) ~ -—V = -1, (5.23) 



n-v y ' V 

1 J? 4 

dx V„5> ~ -. (5.24) 

V — n V R 2 v ; 

Therefore, the criterion of applicability of the late time expansion Q5.10 ) becomes s ^> 
1/V^R 2 . Together with ([513D it yields the new overlap range 

R 2 » a » (5.25) 

which is obviously not empty if the potential satisfies (|5.22|) . 
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To find s* in this case we have to solve the equation (|5.4j ) for the case when s*V is 
not anymore a small quantity. Since V is big, the exponent in ( |5.4j) can be replaced by 
zero inside the compact support, exp(— s*V(x)) ~ 0, |x| < R, and by one outside of it 
where the potential vanishes, exp(— s*V(x)) ~ 1, \x\ > R. Rewriting the integrals in 
both sides of the equation (|5.4| ) as a sum of contributions of \x\ < R and \x\ > R, we 
see that the contribution of the non-compact domain gets cancelled and the equation 
becomes 



J dx. (5.26) 

|a;|<R \x\<R 

Then it follows that is approximately given by the inverse of the function VQ(x) 
averaged over the compact support of the potential 

(5.27) 



/ dxV® 

<v.) s (5 , 8) 

\x\<R 

A qualitative estimate of (V$) ~ Vq implies that s* ~ 1/Vo and it belongs to the 
middle of the interval ( 5.25|) . This makes the case of big potentials fairly consistent. 



On the other hand, the value of $(x) is close to zero inside the potential support 
(see ( 5.23 )), so most likely the estimate for (V$) is smaller by at least one power of 



the quantity 1/VoR 2 , which is the basic dimensionless small parameter in this case. 
Therefore the magnitude of s* becomes bigger by one power of VoR 2 , s* ~ R 2 , which 
is again near the upper boundary of the overlap interval (|5.25 ). Similarly to the small 



potential case, a more rigorous analyses is needed (maybe for more concretely specified 
potentials) to account for subtle edge effects at the boundary of compact support, 
which might shift the value of to a safe region inside ( |5.25|) . 

With the above estimate for ~ R 2 the magnitude of s*V in the expression for 
the infrared part of the effective action fl5.7|) becomes big, s*V ~ s*Vo ~ V R 2 1, 
and we use the big x asymptotics in (|5.8|) to get the contribution 



, d 4 xV$ = —— / dx(VQ) . (5.29) 
64tt 2 s* J 64tt 2 J x 1 v ' 

\x\<R 

In this case the Coleman- Weinberg term is not cancelled in complete agreement with 
what we would expect for big potentials and the final result reads 

r = r div + r cw + ^ / d 4 x(v<s>) 2 . (5.30) 



\x\<R 
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6. Comments 



We developed a new technique for the calculation of late time asymptotics of the heat 
kernel and its functional trace. Using these asymptotics we found previously unknown 
essentially nonlocal and nonperturbative contributions to the effective action for two 
large classes of potentials with compact supports. Therefore, the generalization of 
these results to potentials with power-law falloff, which would imply subtler analyses, 
deserves further studies. 

Our results in their present form are applicable only in higher dimensions d > 3. 
One can easily check that the expression (1/D)V is not well defined in low dimensions 
(d < 2). Note, that the logarithmic kernel of the Green's function of the massless field 
is defined in d — 2 only up to an additive constant. Moreover, the convolution of 
the Green's function with the potential makes sense in d = 2 only if the latter is the 
total derivative of some other function H, as, for instance, the 2-dimensional curvature 
scalar in the Polyakov action M. Thus, the extension of our results to low- dimensional 
models, where other calculational schemes completely fail, is especially important. This 
will be done in the forthcoming paper j7[ . 

Another possibility is the generalization of our technique to potentials with isolated 
zeroes in the interior of spacetime. Even more interesting is the situation when the 
potential becomes negative V(x) < in some spacetime domains. In this case there 
is a tachionic instability, and it is worth deriving quantitative criteria describing this 
instability in terms of the properties of V(x). 

Finally, it is important to generalize our results to curved spacetimes and in addition 
consider fields with non-trivial spin-tensor structure. All these issues are addressed in 
[0. The nonlocal effective action can then be applied to study interesting physical 
problems, like quantum black holes evaporation ||, quantum cosmology, etc. 

It is worth mentioning that the developed technique for late time asymptotics of the 
heat kernel could also be useful in statistical physics for calculating low temperature 
partition functions. 

A. Integral equation for W 

In this Appendix we derive the integral form of the equation 

dW ^ x ' y) + ^yL V w{8\ x, y) - n W {s\ x, y) = f(s\ x, y), (A.l) 

OS s 

where 

f(s\ x, y) = V(x) - (VW(s\ x, y)) 2 , (A.2) 
and W(s — 0| x, y) — 0. With this purpose we first introduce the new function W 

W(s\x,y) = e- sD W(s\x,y). (A.3) 
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Using the relation 

e sn {x - yYe- sn = {x- yY + 2sV^ 
one finds that this function satisfies the following equation 
dW - yY 



ds 



+ 



-V,W = e su f, 



(A.4) 



(A.5) 



which does not contain anymore □-term on the left hand side. To write down the 
formal solution of this equation in terms of the "source" term / = f(s\x,y), let us 
introduce the characteristic curve x^{t) of (|A.5|) , which satisfies the equation 



dx*{t) _ (x(t) - yY 
dt ~~ t ' 

with the boundary conditions 

= 0) = y", x»(t = s)= x". 
The solution of ( |A.6[ ) is 



(A.6) 



(A.7) 



The total derivative of W(t\x(t),y) with respect to t along this characteristic curve is 
then equal to 



jW(t\x(t),y) 



d (x — yY d 
di H t dx^ 



W(t\x,y) = e tn f(t\x(t) : y), 



(A.9) 



where □ = d 2 / dx^dx^. Integrating this equation from to s with the initial condition 
W = at t = and taking into account the boundary conditions ([A .7]) for x(t), one 
gets 



W{s\x,y)= / dte tu f(t\x(t),y) 



(A.10) 



Returning to the original W which is related to W via ( |A.3| ) and taking into account 
that □ = {t/s) 2 □ we finally obtain: 



W{s\x,y) = s / dae sa{1 - a) ^f{sa\x,y) 



x=ax+(l—a)y 



(A.ll) 



where instead of t the new integration variable a = t/s was introduced. This is exactly 
the integral form ( |3.3j ) of the equation ( |A.1|) that we used in Sect. 3. 
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B. Covariant perturbation theory and late time be- 
havior of the functional trace of the heat kernel 



Here we consider the nonlocal covariant perturbation theory (CPT) of [|, |5|. In 
CPT the functional trace of the heat kernel for the covariant second order differential 
operator is expanded as nonlocal series in powers of the potential V with explicitly 
calculable coefficients - nonlocal form factors F n (s\ xi,X2, ■ ■■x n ). Their leading asymp- 
totic behavior at large s was obtained in [|J. Here we calculate them up to the first 
subleading order in 1/s inclusive for a simple operator F — □ — V . Then we explicitly 
perform infinite summation of power series in the potential to obtain the nonlocal and 
nolinear expression (|4.14j) for late time behavior of Tri^(s). 



According to the heat kernel trace is local in the first two orders of the pertur- 
bation theory in the potential ( |2.18| ) 



TrK Q (s) = j—^- 2 Jdx, (B.l) 
TrK 1 {s) = -jjAp 5 j dxV(x), (B.2) 
and in higher orders it reads as 



Tr K n (s) = / J / dx(e sa »)V(x 1 )V(x 2 )...V(x n ] 



, n > 2. 

— X n — X 



(B-3) 



Here Q n is a differential operator acting on the product of n potentials 

71—1 71—1 J— 1 

fl » = E V w + 2 EE &(1 - &) Vi+iV fc+a , (B.4) 

t=l i=2 k=l 

expressed in terms of the partial derivatives labelled by the indices i implying that 
Vj acts on V{xj). It is assumed in ( |B.3| ) that after the action of all derivatives on the 
respective terms all Xi are set equal to x. It is also assumed that the spacetime indices of 
all derivatives V = V M are contracted in their bilinear combinations, VjVfc = VfV^fc. 
The differential operator ( p.4j) depends on the parameters fa, i = 1, ...n — 1, which are 
defined in terms of the parameters a i: i = 1, ...n as 

fa = a i+1 + a i+2 + ... + a n , 

and the angular brackets in (e snn ^ imply that this operator exponent is integrated over 
compact domain in the space of a-parameters 



/71 
d n aS(y^ai - lj exp(sQ n ) 



24 



The late time behavior of Tr K n (s) is thus determined by the asymptotic behavior of 
this integral at s — > oo, which can be calculated using the Laplace method. To apply 
this method, let us note that Q n is a negative semidefmite operator (this is shown 
in the Appendix B of [f|]) which degenerates to zero at n points of the integration 
domain: (0, ...0, cti = 1,0, ...0), i — 1, ...n. Therefore the asymptotic expansion of this 
integral is given by the contribution of the corresponding n maxima of the integrand 
at these points. The integration by parts in (|B.3|) is justified by the formal identity 
Vi + V2 + ••• + V„ = . Using it one can show that the contributions of all these 
maxima are equal, so that it is sufficient to calculate only the contribution of the point 
oil — 1) — 0, i — 2,...n. In the vicinity of this point it is convenient to rewrite 
the expression for Q n in terms of the independent (n — 1) variables 02,^2; ■••«n ) the 
remaining ai — 1 — ^2i=2 a ^ 



^n = ^«iA 2 - a mOikD m D k , (B.5) 

t=2 m,k=2 

where the operator D m is defined as 

D m = V 2 + V 3 + ... + V m , m = 2,.. .n. (B.6) 
Substituting this expression for Q n in (|B.3|) and expanding in powers of the term 



bilinear in a-parameters one gets 

Tr K n {s) = ^j- 2 J dx jy^a exp (s^atf 



1 - s oi m a k D m D k + ... ) ViV 2 ...V n . (B.7) 



m,k=2 



Here 1/n factor disappeared due to the contribution of n equal terms and the range of 
integration over a 2 , ...a n , YH=2 a i — 1> was extended to all positive values of a«. This 
is justified since the error we make by extending the integration range is exponentially 
small and goes beyond the the accuracy of asymptotic expansion in inverse powers of 
s. The second term in the round brackets can be rewritten in terms of the derivatives 
with respect to D 2 m acting on the exponential, so that 



Tr Kr, (s) = — — rrTTT / dx 



( 1 n 

1 - - E D ^ 

V m,k=2 



d d 



U ns )d/2 J I s ^ m K dD 2 m dD 2 h ' 

\ m,k=2 m k 

/ d n ~ l a exp (s ^iD^ V x V 2 ...V n . (B.8) 

•J ■ o 



i=2 



In this form it is obvious that further terms of expansion in powers of the quadratic in 
a part of Q n bring higher order corrections of the 1 / s-series. Doing the integral over 
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a here and performing differentiations one obtains 

f 



TrKJs 



' {4irs) d / 2 

n— 1 n 



+*E E 



da; 
f 



— s 



1 



2 D\...D m _ x (D^) 2 -D^+i- 



D 



m=2 



k ^ +1 Dl-Dl-i (Dl) 2 Dl l+l ...DU (DIY Dl +r ..Dl 



+0[- 

s 



VtV 2 ...V n . 



(B.9) 



The first term in the square brackets gives the leading order term of the late time 
expansion. It can be further transformed by taking into account that any operator D m 
defined by (p.6| ) acts as a partial derivative only on the group of factors V m V m+ i...V n 
in the full product Vi...V n , D m Vi...V n = Vi...V m -i S/(V m +i...V n ). Therefore all the 
operators understood as acting to the right can be ordered in such a way 



TrKJs) 



( 47rs )rf/2 



dxV 1 -^V n ^V n _ 1 ...-^V 2 + ' ' 



D 2 D 2 

u n u n-\ 



d/2 



that the labels of -D^'s can be omitted and all can be identified with boxes also 
acting to the right 



TrKJs) 



( 4vrs )rf/2 



dx V-V-...V- V(x) + ' 
□ □ □ v ' \s 



d/2 



(B.10) 



n-l 



Infinite summation of this series is not difficult to perform because this is the geometric 
progression in powers of the nonlocal operator V(l/n) and 



or 



TiK(s) = TrK (s) 



TtK(s) 



(4vrs) d / 2 



dxil-sD 



1 



□ - V 



v(x)+o(s°: 



(B.11) 



{A7Ts) d / 2 

The second term here looks as a total derivative. However, it does not vanish because 
this is a derivative of the nonlocal expression and the corresponding surface term does 
not vanish at infinity in view of the Green's function asymptotics. This term can be 
rewritten as 

1 „, , .„ , „ 1 



□ - 



□ - V 
and, therefore, 



V(x) = V(x) + V 



n-v 



V(x) = V$(x) 



TrK(s) = Jdx(l- sV$(x) + 0( a °)) : 



(B.12) 



(B.13) 
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where 0(s°) denotes the subleading in s terms which depend on the potential in a 
nontrivial way. They are given by infinite resummation over n of the second and 
third terms in square brackets of Eq. ( P-9|) . Remarkably, this summation can again be 
explicitly done. In this case one has to sum multiple geometric progressions. 
Indeed, the second term of (|B.9|) gives rise to the series 



( 4vrs )d/2 



co n 111 1 

dxYYv -V...-V -z V-...V- V(x). (B.14) 

J ^ — O ™ — o \ ✓ \ ✓ 



n=2 m=2 

n—m m—2 



By summing the two geometric progressions with respect to independent summation 
indices < n — m < oo and < m — 2 < oo one finds that this series reduces to 

dxV- ——V(x), (B.15) 



(4vrs) rf / 2 J (□ - V) 2 

which after the integration by parts amounts to 



2 



Similarly, the third term of ( |B.9|) gives rise to the triplicate geometric progression 
which after summation and integration by parts reduces to 



^ \ n J a a ^ \ a J □ □ ^ \ □ 

i=0 j=0 z=o 

1 . .1 



(v M $(x)) ——V-V^(x). (B.17) 



(47rs) d / 2 7 V M v 7 □ - K □ 

Taking here into account that 

1111 

V- = (B.18) 



(47rs) d / ; 



□ - V □ □ - V □ 
one finds that the sum of (|B.16| ) and (|B. 17] ) is equal to 

where the cancellation of the first and the third terms takes place after rewriting 
in the third term as V M ($ — 1) and integrating it by parts[]. Together with ( P-13|) the 
contribution ( |B.19 ) forms the nonlinear and nonlocal late time expression for the heat 
kernel trace ( f4.14|) up to the first subleading order in 1/s inclusive. 

9 Straightforward integration by parts of V M <f>(l/D)V M $ is impossible because $(x) does not vanish 
at |x| — > oo, while $(x) — 1 docs. 
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